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Abstract 

To find all extreme points of multimodal functions is called extremum problem, which is a well 
known difficult issue in optimization fields. Applying ant colony optimization (ACO) to solve this 
problem is rarely reported. The method of applying ACO to solve extremum problem is explored 
in this paper. Experiment shows that the solution error of the method presented in this paper is 
less than 1CP 8 . 
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I. INTRODUCTION 



Extremum Problem 



Multimodal function refers to the function which has more than one extreme point. 
To find all extreme points is called extremum problem, which is a well known difficult 
issue in optimization fields. Many practical engineering problems can be converted as this 
problem, such as the detection of multiple objects in military field. Therefore, solving the 
extremum problem is a useful study topic. To solve the extremum problem, many methods 
of optimization are applied, such as genetic algorithm (GA) simulated annealing (SA) 
[2I, particle swarm optimization algorithm (PSO) js, 4], immune algorithm ( IA) and so 
on. However, currently there is rare report that applying Ant Colony Optimization (ACO) 
to solve the extremum problem. The motivation of this paper is to apply ACO to search all 
extreme points of function. 



B. Introduction of Ant Colony Optimization (ACO) 

Ant Colony Optimization (ACO) was first proposed by Dorigo (1991) The 
inspiring source of ACO is the foraging behavior of real ants. When ants search for food, 
they initially explore the area surrounding their nest in a random manner. As soon as an 
ant finds a food source, it remembers the route passed by and carries some food back to 
the nest. During the return trip, the ant deposits pheromone on the ground. The deposited 
pheromone, guides other ants to the food source. And the feature has been shown, indirect 
communication among ants via pheromone trails enables them to find the shortest routes 
between their nest and food sources. ACO imitates this feature and it becomes an effective 



algorithm for the optimization problems 



1 



It has been successfully applied to many com- 



binatorial optimization problems 10, llll, ll2|], such as Traveling Salesman Problem (TSP) 
13, Ijj], Quadratic Assignment Problem(QAP) If] ]. Job-shop Scheduling Problem(JSP) 



16), Vehicle Routing Problem(VRP) (l7|,ll8j, Data Mining(DM) [19] and so on. 



The application of ACO pushes the study of ACO theory, and its two main study topics 
are the analysis of convergence and runtime. M. Birattari proves the invariance of ACO 
and introduced three new ACO algorithms 20]. Convergence is one of focus study of ACO. 



Walter J. Gutjahr studied the convergence of ACO firstly in 2000 



2l|. T. Stwezle and 



M. Dorigo proved the existence of the ACO convergence under two conditions, one is to 
only update the pheromone of the shortest route generated at each iteration step, the other 



is that the pheromone on all routes has lower bound 22j. C.-Y. Pang and et.al. found a 
potential new view point to study ACQ convergence under general condition, that is entropy 
convergence m 2009 [J. Inref.lJ, the following conclusion is get: ACO may not converges 
to the optimal solution in practice, but its entropy is convergent. The other study focus 
of ACO is time complexity. ACO has runtime 0(tmN 2 ), where t, m and iV refers to the 
number of iteration steps, ants, cities and M = [j^] in general. To reduce runtime, cutting 



down parameter t and N is the main way possi 



some theoretical results about ACO runtime 
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y. In 2008, Walter J. Gutjahr presented 



251 ]. Since runtime is proportional to the 



square of N, parameter N is the key factor of runtime. Through cutting down N to reduce 
runtime is a choice. And Pang and et al. do the following attempt [26]. Firstly, cluster all 
cities into some classes (group) and let ACO act on these small classes respectively to get 
some local TSP routes. And then joint these local route to form the whole TSP route. If 
class is compact, the length of local route got at every iteration step will change continually 
possibly, where compactness refers to all data cluster in a small region tightly. And this 
property result in the conclusion: the convergence criterion \ Lt ~Lt+i\ q is the marker 
of ACO convergence, where L t is the length of local route at t — th iteration step. Thus, 
minimum iteration number t can be estimated by the marker approximately. 

The study of ACO theory speeds up its application again. ACO not only can be applied 
to solve discrete optimization problems, but also to continuous ones. The first method for 
continuous-space optimization problems, called Continuous ACO (CACO), was prop osed by 
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33|. In 



Bilchev and Parmee (1995) [29], and later it was used by some others 
general, the application of ACO to continuous optimization problems need to transform a 
continuous search space to a discrete one. Other methods include that, in 2002, Contin- 
uous Interacting Ant Colony (CIAC) was proposed by Dreo and Siarry 34j, and in 2003, 
an adaptive ant colony system algorithm for continuous-space optimization problems was 



proposed by Li Yan-jun 



and so on. 
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C. Framework of ACO 



Traveling Salesman Problem (TSP) is a famous combinatorial problem, it can be stated 
very simply: 

A salesman visit N cities cyclically provided that he visits each city just once. In what 

order should he visit them to minimis the distance traveled? 

The typical application of ACO is to solve TSP, and its basic idea is stated as below: 
When an ant passes through an edge, it releases pheromone on this edge. The shorter 

the edge is, the bigger the amount of pheromone is. And the pheromone induces other ants 

to passes through this edge. At last, all ants select a unique route, and this route is shortest 

possibly. 

The framework of ACO is introduced as below: 

Stepl (Initialization): Posit M ants at different M cities randomly; Pre-assign maxi- 
mum iteration number t max ; Let t — 0, where t denotes the t — th iteration step; Initialize 
amount of pheromone of every edge. 



Step2.1: Every ant selects its next city according to the transition probability. The 
transition probability from the i — th city to the j — th city for the k — th ant is defined as 



, where allowed], denotes the set of cities that can be accessed by the k — th ant; ^(t) is 
the pheromone value of the edge is local heuristic function defined as 

1 

, where is the distance between the i — th city and the j — th city; the parameters a 
and (3 determine the relative influence of the trail strength and the heuristic information 
respectively. 



Step2.2: After all ants finish their travels, all pheromone values are updated according 
to EqfJJ 



Step2 While(t < t max ) 
{ 



Eqffl 




(1) 
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T ij (t+l) = (l- P )-T ij (t)+AT ij (t) (2) 
m 

Ar y (*) = X)ArW(t) 
fc=i 

if the k — th ant pass edge 
0, otherwise 

, where (t) is the length of the route passed by the k — th ant during the t — th iteration; 
p is the persistence percentage of the trail (thus, 1 — p corresponds to the evaporation); Q 
denotes constant quantity of pheromone. 

Step2.3: Increase iteration step: t <— t + 1 

} 

Step3: End procedure and select the shortest route as output from the routes traveled 
by the ants. 

II. APPLY ACO TO SEARCH ALL EXTREME POINTS OF FUNCTION 
A. Basic Idea 

Assume that the function is f(x), x is real number and it belongs to a closed interval 
[a, b]. The task of this paper is to extract all extreme points that the corresponding value is 
minimal locally The basic idea of this paper is stated roughly as below: 

Divide interval [a, b] into many tiny intervals with equal size. Suppose these small interval 
are {Ii,I 2 , •■• ,I n } and the center of interval ij is denoted by Xj. Suppose the neighbor 
interval of interval ij is ij+i (or U-i). And an ant is put at the center of each small interval. 
If f(xi) > f(x i+ i), the ant at interval ij will move to interval possibly, just liking there 
is an virtual edge between Jj and And assume the virtual edge is e(ij,/ i+ i). The 

weight (virtual distance) of edge e(Ii,I i+ i) is proportional to f{xj) — f(x i+ i). That is, the 
bigger f(xi) — /(xj+i) is, the more possibly the ant moves to from Jj. When the ant 
moves to Ii+i, it releases pheromone at ij+i, and the pheromone is proportional to value 
f{ x i) ~ f( x i+i)- The pheromone depositing on will attract other ants move to it. 
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After some iteration steps, some intervals contain many ants while other ones contain no 
ant. The intervals containing ants include extreme points possibly, while other ones include 
no extreme point possibly. And then keep the intervals which contain ants, and divide them 
into much more small intervals, repeat the same procedure again until the size of intervals 
is sufficient small. At last all ants will stay around extreme points. The centers of these 
sufficient small intervals are the approximations of extreme points. 

From above discussion, it can be seen that the realization of the basic idea consists of four 
parts: partition of interval [a, b] and initialization, rule of ant moving, rule of pheromone 
updating, and keeping the intervals containing ants. The contents of the four parts are 
stated as below. 

B. Partition of Interval and Initialization 

Suppose interval [a, b] is partitioned into n small intervals with equal size, which are 
denoted by {I±, I 2 , • • ■ , I n }, where n is a pre-assigned number. Then each interval has size 

b — a 



Ii = [a + {i — 1)5, a + iS] 
The i — th interval has center Zj. 

Xi = a + {i- -)6 

Suppose t denotes the t — th iteration step of ACO and it is initialized as zero (i.e., 
t = 0). Put n ants at the centers of the n intervals, and each interval has only one ant. 
Suppose these ants are denoted by 01,02,- • • a n respectively, and ant is associated with 
the i — th interval Jj. Each ant will release an initial pheromone at its associated interval 
Ii (i.e., Tj(0) = const, const is a constant number). In addition, set the increment of the 
pheromone of each interval to zero (i.e., Arj(0) = 0) 

Figure 1 shows a diagram of initialization. 
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FIG. 1: Initialization: Divide the function domain into intervals and put an ant at the center of 
every interval. 

C. Rule of Ant Moving 

Let Neighbor (Ii) be a set of neighboring intervals of Ij. Take one-dimensional function 
for example 



Neighbor(Ii) 



Ui+i} , i = 1 

{/<_!, J i+ i}, i = 2---n-l 
{Ii-i} , i = n 

As it is discussed in section A, the ant staying at interval Ii will move to neighbor interval 
denoted by Ij, just like there is a virtual edge e(Ii,Ij). The weight of the virtual edge is 
1/(^0 — f( x j)\- Then the heuristic factor is 

Vij = \f( x i) ~ f( x j)\ 

Suppose interval Ii contains ant a&. If > f(xj), ant is allowed to move to neighbor 
interval Ij. Otherwise, it is forbidden to move. Suppose all intervals allowed to be accessed 
by ant is marked as allowedk- The transition probability of ant cik is defined as 

T ' [! '" i - -jr, if f{ Xi ) > f( Xj ) 



Pij (.*) = < (Oj 

0, otherwise 
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In Eq|3j a is the relative influence of the trail strength; (3 is the heuristic information; 
Tj(t) is the pheromone of interval Ij. 

D. The Rule of Pheromone Updating 

Suppose dk ant is staying at interval li and it will move to neighbor interval Ij. After it 
moves to Ij, releases pheromone at Ij. The pheromone amount is denoted by Arj k \t) 

Arj h \t) = C 1 (f{x i )-f{x j )) 
, where G\ is a positive constant. 

The bigger f(xt) — f(xj) is, the higher the amount of released pheromone is, the more 
possibly that other ants will be attracted to interval Ij. 

Not only ant at arrives Ij and releases pheromone, but also other ants which move to 
interval Ij and release pheromone too. Suppose there are q ants will move to interval Ij 
during the t — th iteration step, which are denoted by The sum of released 

pheromone by aj 1 , a j2 , • ■ ■ , dj q is At, (t) . Then 

A Ti (i) = X>^(t) 
P =i 

When all ants a^, % 2 , • * * > a j q m ove to Ij, the amount of pheromone Tj is changed as 
EqJl 

T j (t+l) = (l-p)-T j (t)+A Tj (t) (4) 
, where p is the evaporation percentage of the trail (thus, 1— p corresponds to persistence). 

E. Keeping Only the Intervals Containing Ants to Cut Down Search Range 

The intervals that have smaller function values depositing much more pheromone, and 
it will attract ants more powerfully. After several iterations, the distribution of ants has 
the feature that all ants stay at the intervals that have smaller function values and other 
intervals contain no ants. That is, extreme points are included in the intervals containing 
ants. And then keep the intervals having ant and delete other intervals to update search 
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range. Thus, the updated search range became smaller. Dividing the updated search range 
into smaller intervals will result in much smaller search range at next iteration step. When 
intervals becoming sufficient small, all ants will stay around extreme points, the centers of 
the intervals are their approximations. 

Figure 2 shows the distribution feature of ants. 

In addition, as it is well known, ACO runs slow, which is the bottleneck of application. 
And it evades this bottleneck that keeping only the intervals having ants to cut down search 
range. 




FIG. 2: The Feature of Ant Distribution: The intervals including extreme points have more 
pheromone to attract ants. After several iteration steps, all ants will stay at the intervals con- 
taining extreme points while other intervals are empty. Then search range become smaller, which 
consists of the intervals containing ants. At last all ants stay around extreme points. 
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F. Method of Searching All ExtremePoints 

Stepl (Initialization): Divide domain [a, b] into many small intervals and put an ant 
in each interval; Do other initialization. The detail is shown at section B. Suppose 5 is the 
length of interval and e is a stop threshold. 

Ste P 2 While (5 > e) 

{ 

Step2.1: All ants move to new intervals according the rule shown at section C 
Step2.2: Update pheromone according the rule shown at section D 
Step2.3: Update search range according to section E and divide it into smaller intervals 
(Suppose the number of these intervals is nl) . Calculate the size of interval and set it to S. 

} 

Step3 Extract all intervals that contain ants, the centers of the intervals are the approx- 
imations of extreme points. 

If argument x is multi-dimensional vector, divide the range of every component of vector 
into smaller intervals, the combination of these intervals forms many small lattices. And 
then put an ant in each lattice, apply the above method, all extreme points can be extracted. 

G. An Example 

To understand above method easily, an simple example is stated as below: 

Assume that the domain of 1-dimensional function is divided into 3 intervals I\, J 2 , I3, 
which associated center is Xi, x 2 and £3 respectively. Initially ant a±, a 2 , and 0,3 is put at 
xi, X2 and 23 respectively. 

Check the first ant: If f(x\) > /(x 2 ), ant a\ moves to interval J 2 . Otherwise, do nothing. 

Check the 2nd ant: If their is unique interval (e.g. I3 ) such that /(a^) > f{ x 3), ant a 2 
moves to interval I3. If /(x 2 ) is smaller than both f(x\) and f{x^), do nothing. If /(a^) is 
bigger than both f{x\) and f(xs), it is uncertain that ant a 2 moves to I\ or J 3 . And ant 
a 2 will select its visiting interval randomly according to its transition probability defined at 
EqU 

Check the 3rd ant using same way. 

After all ants are checked, update their associated interval (position) and interval 
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pheromone. Repeat above processing until all ants can not move. 

Then keep the intervals which contains ants, and delete other blank intervals. And 
divided the intervals containing ants into smaller interval, repeat above process until the 
size of interval is sufficient small. And then all interval centers are the approximations of 
extreme points. 

III. EXPERIMENT 

In this section, several functions will be tested. The parameters are listed as below: 
const = 10, a = 1, f3 = 1, C x = 1, p = 0.3, e = 0.0001 

Two performances are considered, which are error (ratio of inaccuracy) r and runtime. 
Error r is defined as 



, where (xo,f(xo)) denotes the true extreme point on theory and (x' ,f(x' )) is its ap- 
proximation calculated by the method presented in this paper. 

In addition, the hardware condition is: notebook PC DELL D520, CPU 1.66 GHZ. 

A. Instance 1 (see table 1 and Figj3]): 

fi(x) = sin 6 (5.l7rx + 0.5), x G [0, 1] 

In the experiment, additional parameter is n = 20 and nl = 10. Table 1 and Figj3]show 
all extreme points (local maximal points) of fi{x). Only 0.7106 seconds is cost. 

Other functions are tested, their errors are less than 10 -8 except the boundary. And 
runtime is less than 1 second (see appendix I). 

B. Instance 2 (see Figj4j): 



Instance 2 is a typical test function, which include many extreme points and any small 
change of argument x will result in big change. In addition, the theoretical calculation of 
extreme points of instance 2 is difficult. 




f 2 (x) = 5e 



-0.5x 



sin(30x) + e°- 2x sin(20x) +6, x G [0, 8] 
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The additional parameters are n = 480 and nl = 10. Fig|4] shows all the calculated 
extreme points, and the real numbers are listed at appendix (see appendix II). 

C. Instance 3 (see Figj5]): 

f3(x±, X2) = x\ + x\ — cos(18xi) — cos(18x2), X\,X 2 G [—1, 1] 

The interval [— 1, 1] ® [— 1, 1] is divided into n = 40 x 40 intervals initially. And at next 
iteration steps, search domain is divided into nl = 20 x 20 small intervals. 36 extreme points 
are got and shown at FigJSl and the digital solutions are listed at appendix (see appendix 
III) 

Instance 3 is a 2-dimensional function and 3203.2968 seconds is cost. And it is slower 
than 1-dimensional function instance 1 and instance 2. To improve running speed is the 
next work. 

Many functions are tested by the authors, and some experiments results are listed at 
appendix. These tests demonstrate that solution error is less than 10~ 8 except the special 
case that extreme point is at the border of the domain. In addition, these testes also 
demonstrate that the method is very fast for 1-dimensional function. 
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TABLE 1. The Extreme Points of Function fi(x) 


Extreme 


Theory Value 


Calculated Value 


Error (%) 


\- i "(l Tl t" ^ 


1 Tn / 1 Tr, 1 1 

^0) J K^o)) 


(V f(r')) 


1 


1 


0.066832364, 
i 

± 


0.066827175, 
n qqqqqqq7q 


2.1e-06 


2 


0.262910795, 
i 

± 


0.262914975, 
n qqqqqqqs7 


1.3e-06 


3 


0.458989227, 
i 

± 


0.458990475, 
n qqqqqqqqq 


1.0e-07 


4 


0.655067658, 
i 

± 


0.655066175, 
n qqqqqqqqs 


2.0e-07 


5 


0.851146090, 
1 


0.851148550, 
0.999999995 


5.0e-07 


6 


1 

0.147822118 


0.999998400, 
0.147800654 


0.0145 


Conclusion: Error of the method presented this paper is low except the border point ( 6th point). 
And runtime is fast (0.7106s). Notice: The 6 points are shown at Figj3] 



Notice: From the Table 1, we can see that the border point has big error because the 
value calculated is the center of the interval, not boundary. To evade this drawback, the 
function value at boundary can be calculated directly. 



IV. CONCLUSION 

To find all extreme points of multimodal functions is called extremum problem, which 
is a well known difficult issue in optimization fields. It is reported rarely that applying 
Ant Colony Optimization(ACO) to solve the problem. And the motivation of this paper 
is to explore ACO application method to solve it. In this paper, the following method is 
presented: 

Divide the domain of function into many intervals and put an ant in each interval. And 
then design rule such that every ant moves to the interval containing extreme point near by. 
At last all ants stay around extreme points. 

The method presented in this paper has following three advantages: 
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FIG. 3: The Extreme Points Calculated of Function fi(x). 




FIG. 4: The Extreme Points Calculated of Function f2(x) 

1. Solution accuracy is high. Experiment shows that solution error is less than 10 -8 . 

2. Solution calculated is stable (robust). Ant only indicates the interval containing 
extreme point, not the accurate position of extreme point. It is easy for ant to find a 
interval although finding a special point in interval is difficult. 

3. The method is fast for 1-dimensional function. ACO is slow. But some feature is 
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FIG. 5: The Extreme Points Calculated of Function f3(xi,xz). 
found to speed ACO (see section 2.5) 
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V. APPENDIX 
A. Appendix I 
1. Instance 4- 

U(x) = (x + l)(x + 2)(x + 3)(x + 4)(x + 5) + 5 x E [-5, 0] 
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The additional parameters are n = 30 and ri\ = 10. 



Table 2. Extreme Points of Function fi(x) 










Extreme 


Theory Value 


Calculated Value 


Error (%) 


Points 


(x ,f(x )) 




r 


1 


-5.00000000000000, 
5.00000000000000 


-4.99999893333333, 
5.00002559994310 


5.12e-004 


2 


-3.54391225590234, 
3.58130337441708 


-3.54391549166667, 
3.58130337448565 


1.91e-009 


3 


-1.35556713184173, 
1.36856779155116 


-1.35557000833333, 
1.36856779171500 


1.20e-008 


Notice: 1. Runtime is 0.5280s. 2. From the table, we can see that only the boundary 


maximum value has big error because x' is the center of interval, not the boundary. 


To evade this drawback, the function value at boundary can be calculated directly. 





140 




FIG. 6: The Extreme Points Calculated of Function fi(x). The red stars are the extreme 
points. 
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2. Instance 5: 



f 5 (x) = (x + 2) cos(9x) + sin(7x) x £ [0, 4] 

The additional parameters are n = 95 and n\ = 10 



Table 3. Extreme Points of Function fs(x) 










Extreme 


Theory Value 


Calculated Value 


Error (%) 


Points 


(x , f(x )) 


\ U i v \ U/ / 


r 


1 


0.0000000000000, 
2.0000000000000 


0.00000134736842, 
2.00001077880032 


5.39e-004 


2 


0.38762423574739, 
-1.8300394960224 


0.38762006315789, 
-1.83003949456305 


7.97e-008 


3 


1.03491254831578, 
-2.19650045294125 


-1.03490865263158, 
-2.19650045140718 


16.98e-008 


4 


1.72787330641709, 
-4.13597012121530 


1.72787490526316, 
-4.13597012080927 


9.82e-009 


5 


2.44888001781347, 
-5.43427465397202 


2.44888418947368, 
-5.43427465041009 


6.55e-008 


6 


3.16064032530539, 
-5.21793460286144 


3.16064037894737, 
-5.21793460286083 


1.17e-011 


7 


3.84493263713282, 
-4.86068863199737 


3.84493433684211, 
-4.86068863138223 


1.27e-008 


Notice: 1. The runtime is 0.9030s. 2. The boundary maximum value has big error 


because x' Q is the center of the interval, not the boundary. To evade this drawback, 


the function value at boundary can be calculated directly. 





B. Appendix II. 

f 2 (x) = 5e-°- 5x sin(30a;) + e°- 2x sin(20x) + 6, x G [0, 8] 
The additional parameters are n = 480 and nl = 10. 
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Table 4. Extreme Points of Function f2(x) 




n ."vi" TPTTl P 
J_JjV LI C1IIC 


CXiLKj LI Id LCVJ. V C11U.C 


T-rvf tptti p 
J_JjV LI C1IIC 


( , n 1 ni 1 pi "hf^ri x/aliif^ 

CXlLy LllcL LCLL V CX1U.C 


n "vi" ypm p 

J_JjV LI CIIIC 


I , PI 1 f*11 1 5l1"pH 

V_y CllL/ LLlCL LCwl V CXIU.C 


Points 


K,/K)) 


Points 


(x' Q , f(x Q )) 


Points 


(x , f(x Q )) 


1 


n nnnn 

6.0005 


12 


5.5878 


23 


5 9547 
3.3639 


2 


U. ±U-LtJ, 

1.3372 


13 


9 4^91 
3.0483 


24 


^ ^819 
2.7148 


3 


2.7441 


14 


9 (SQ89 
4.2229 


25 


5 88(S9 
2.9274 


4 


1.2573 


15 


9 83IS8 
5.9928 


26 


fi 9081 
2.3729 


5 


2.5460 


16 


3 0777 
3.2186 


27 


<S ^1 (S4 
2.4487 


6 


n Q88Q 
3.9745 


17 


o .oouo , 

4.2595 


28 


IS 83^ 
1.9567 


7 


2.1270 


18 


3 7039 
3.2478 


29 


7 1 48^. 
1.9215 


8 


1 49^9 
3.3895 


19 


3 Q899 
4.0777 


30 


7 4<S31 

1 .4:1)01 , 

1.4650 


9 


1 (SI ^ 
4.9025 


20 


4 3988 
3.1621 


31 


7 7748 
1.3366 


10 


1.8261, 
2.7009 


21 


4.6206, 
3.7560 


32 


8.0000, 
7.1737 


11 


2.0601, 
3.9371 


22 


4.9548, 
2.9802 






Notice: Runtime is 4.3411s. The associated figure is FigJJ] 



C. Appendix III. 

fs{x\, X2) = x\ + x\ — cos(18xi) — cos(18x2), Xi, X2 G [—1, 1] 

The interval [— 1, 1] <S> [—1, 1] is divided into n = 40 x 40 intervals initially. And at next 
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iteration steps, search domain is divided into nl = 20 x 20 small intervals. 



Table 5. The Extreme Points of Function fz(xi,x 2 ) 












Extreme 


Calculated Value 


Extreme 


Calculated Value 


Extreme 


Calculated Value 


Points 


(xi,x 2 ,f(xxX 2 )) 


Points 


(x 1 ,x 2 , f{xix 2 )) 


Points 


(x 1 ,x 2 ,f(xix 2 )) 


1 


-0.8781,-0.8781, 
3.5326 


13 


-0.1810, -0.8725, 
2.7873 


25 


0.5175, -0.8700, 
3.0176 


2 


-0.8737,-0.5310, 
3.0362 


14 


-0.1773, -0.5252, 
2.3056 


26 


0.5200, -0.5200, 
2.5367 


3 


-0.8725,-0.1810, 
2.7873 


15 


-0.1756, -0.1756, 
2.0613 


27 


0.5215, -0.1700, 
2.2969 


4 


-0.8700,0.1675, 
2.7759 


16 


-0.1715,0.1715, 
2.0559 


28 


0.5252,0.1773, 
2.3056 


5 


-0.8700,0.5175, 
3.0176 


17 


-0.1700,0.5215, 
2.2969 


29 


0.5268,0.5268, 
2.5517 


6 


-0.8675,0.8675, 
3.4966 


18 


-0.1675,0.8700, 
2.7759 


30 


0.5310,0.8737, 
3.0362 


7 


-0.5310,-0.8737, 
3.0362 


19 


0.1675, -0.8700, 
2.7759 


31 


0.8675, -0.8675, 
3.4966 


8 


-0.5268,-0.5268, 
2.5517 


20 


0.1700, -0.5215, 
2.2969 


32 


0.8700, -0.5175, 
3.0176 


9 


-0.5252,-0.1773, 
2.3056 


21 


0.1715, -0.1715, 
2.0559 


33 


0.8700, -0.1675, 
2.7759 


10 


-0.5215,0.1700, 
2.2969 


22 


0.1756,0.1756, 
2.0613 


34 


0.8725,0.1810, 
2.7873 


11 


-0.5200,0.5200, 
2.5367 


23 


0.1773,0.5252, 
2.3056 


35 


0.8737,0.5310, 
3.0362 


12 


-0.5175,0.8700, 
3.0176 


24 


0.1810,0.8725, 
2.7873 


36 


0.8781,0.8781, 
3.5326 


Notice: All points are shown at Figj5j Runtime is 3203.2968s. 
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